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Abstract 

We study the multifractal temporal scaling properties of river discharge and precip- 
itation records. We compare the results for the multifractal detrended fluctuation 
analysis method with the results for the wavelet transform modulus maxima tech- 
nique and obtain agreement within the error margins. In contrast to previous stud- 
ies, we find non-universal behaviour: On long time scales, above a crossover time 
scale of several months, the runoff records are described by fluctuation exponents 
varying from river to river in a wide range. Similar variations are observed for the 
precipitation records which exhibit weaker, but still significant multifractality. For 
all runoff records the type of multifractality is consistent with a modified version of 
the binomial multifractal model, while several precipitation records seem to require 
different models. 



The analysis of river flows has a long history. Already more than half a cen- 
tury ago the engineer H. E. Hurst found that runoff records from various 
rivers exhibit 'long-range statistical dependencies' [1]. Later, such long-term 
correlated fluctuation behaviour has also been reported for many other geo- 
physical records including precipitation data [2,3], see also [4]. These original 
approaches exclusively focused on the absolute values or the variances of the 
full distribution of the fluctuations, which can be regarded as the first mo- 
ment F\(s) [1-3] and the second moment F 2 (s) [5], respectively. In the last 



Preprint submitted to Elsevier Science 



April 5, 2003 



15000 




1873 1874 1875 
time [years] 



1873 1874 1875 
time [years] 



Fig. 1. Three years of (a) the daily runoff record of the river Danube (Orsova, 
Romania) and (b) of the daily precipitation recorded in Vienna (Austria) . 



decade it has been realized that a multifractal description is required for a full 
characterization of the runoff records [6,7]. Accordingly, one has to consider 
all moments F q (s) to fully characterize the records. This multifractal descrip- 
tion of the records can be regarded as a 'fingerprint' for each station or river, 
which, among other things, can serve as an efficient non-trivial test bed for 
the state-of-the-art precipitation-runoff models. 

Since a multifractal analysis is not an easy task, especially if the data are 
affected by trends or other non-stationarities, e.g. due to a modification of 
the river bed by construction work or due to changing climate, it is useful to 
compare the results for different methods. We have studied the multifractality 
by using the multifractal detrended fluctuation analysis (MF-DFA) method 
[8] (see also [9,10]) and the well established wavelet transform modulus max- 
ima (WTMM) technique [11,12] and find that both methods yield equivalent 
results. Both approaches differ from the multifractal approach introduced into 
hydrology by Lovejoy and Schertzer [6,7]. 

We analyze long daily runoff records {Wi} from six international hydrological 
stations and long daily precipitation records {Pi} from six international me- 
teorological stations. The stations are representative for different rivers and 
different climate zones, as we showed in larger separate studies [13,14]. As 
a representative example, Fig. 1 shows three years of the runoff record of 
the river Danube (a) and of the precipitation recorded in Vienna (b). It can 
be seen that the precipitation record appears more random than the runoff 
record. To eliminate the periodic seasonal trend, we concentrate on the depar- 
tures 4>i = Wi — Wi (and fa = Pi — Pi) from the mean daily runoff Wi. Wi 
is calculated for each calendar date i, e.g. 1st of April, by averaging over all 
years in the record. 

In the MF-DFA procedure [8], the moments F q (s) are calculated by (i) in- 
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Fig. 2. Comparison of the fluctuation functions F q (s) calculated with the multi- 
fractal detrended fluctuation analysis (MF-DFA, filled symbols) with the rescaled 
wavelet transform modulus maxima (WTMM) partition sums [sZ(q,s/8)] 1 / q (open 
symbols) as function of time scale s (in days) for (a) the river Weser (Vlotho, 
Germany, 171y) and (b) the river Danube (Orsova, Romania, 151y). The different 
symbols indicate different moments, q = —6 (triangles up), q = —2 (circles), q = 2 
(squares), q = 6 (triangles down), and the curves are shifted vertically for clarity 
The slopes h(q) for large s of both, the MF-DFA curves and the rescaled WTMM 
curves are equivalent. 

tegrating the series, (ii) splitting the series into segments of length s, (iii) 
calculating the mean-square deviations F 2 (u,s) from polynomial fits in each 
segment, (iv) averaging [F 2 (v, s)] 9//2 over all segments, and (v) taking the gth 
root. In the paper, we have used third order polynomials in the fitting pro- 
cedure of step (iii) (MF-DFA3), this way eliminating quadratic trends in the 
data. We consider both, positive and negative moments F q (s) (q ranges from 
—10 to +10) and determine them for time scales s between s = 5 and s = N/5, 
where is the length of the series. Figure 2 shows the results (filled symbols) 
for two representative hydrological stations. On large time scales, above a 
crossover occurring around 30-200 days, we observe a power-law scaling be- 
haviour, 



F g (s) ~ s h ^\ 



where the scaling exponent h(q) (the slope in Fig. 2) explicitly depends on the 
value of q. This behaviour represents the presence of multifractality. 

In order to test the MF-DFA approach we have applied the well-established 
WTMM technique, which is also detrending but based on wavelet analysis 
instead of polynomial fitting procedures. For a full description of the method, 
we refer to [11,12]. First, the wavelet-transform T(n,s') = 4 Y,iLi 4>i 9[{i ~ 
n)/s'} of the departures 0j is calculated. For the wavelet g{x) we choose the 
third derivative of a Gaussian here, g(x) = d 3 (e~ x ^ 2 )/dx 3 , which is orthogonal 
to quadratic trends. Now, for a given scale s', one determines the positions 
rii of the local maxima of \T(n,s')\, so that \T{n,i — l,s')\ < \T(rii,s')\ > 
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Fig. 3. The generalized Hurst exponents h(q) for six representative daily runoff 
records: (a) Amper in Fiirstenfeldbruck, Germany, (b) Weser in Vlotho, Germany, 
(c) Susquehanna in Harrisburg, USA, (d) Wertach in Biessenhofen, Germany, (e) 
Danube in Orsova, Romania, and (f) Niger in Koulikoro, Mali. The h(q) values have 
been determined by straight line fits of F q (s) on large time scales. The error bars of 
the fits correspond to the size of the symbols. The lines are obtained by fits of the 
two-parameter binomial model yielding Eq. (3). The resulting model parameters a 
and b are reported in the figures. All fits are consistent with the data within the 
error bars (from [13]). 

\T(rii + l, s')\. Then, one obtains the WTMM partition sum Z q (s') by averaging 
\T(rii,s')\ q for all maxima rij. An additional supremum procedure has to be 
used in the WTMM method in order to keep the dependence of Z(q, s') on s' 
monotonous [12]. The expected scaling behaviour is Z(q,s') ~ (s') T ( g \ where 
r(q) are the Renyi exponents. Since r(q) is related to the exponents h(q) by 
h(q) = [r(q) + l\/q [8] we have plotted 
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We set s' = s/8 in the comparison with the MF-DFA results, since the wavelet 
we employ can be well approximated within a window of size 8s' (i.e. within 
4 standard deviations on both sides), and this window size corresponds to the 
segment length s in the MF-DFA. Figure 2 shows that both methods yield 
equivalent results for the q values we considered. 

Using the MF-DFA results, we have determined h(q) from Eq. (1) for all 
runoff records and all precipitation records and for several values of q. Since 
a crossover occurs in F q (s) for time scales in the range of 30-200 days, we 
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Fig. 4. The generalized Hurst exponents h(q) for six representative daily precipita- 
tion records: (a) Arhangelsk (Russia), (b) Hamburg (Germany), (c) Winnemucca 
(USA), (d) Cheyenne (USA), (e) Vienna (Austria), (f) Moskow (Russia), analogous 
with Fig. 3. While the fits in (a,b,d,e) are consistent with the data within the error 
bars, significant deviations occur in (c) and - even more drastically - in (f). 

considered only sufficiently long time scales (above one year), where the results 
scale well. Figure 3 shows h(q) for the runoff data, while Fig. 4 shows h(q) 
for the precipitation data. Together with the results we show least-square fits 
according to the formula 



%) = 



1 In [a 9 + b q ] 
q q In 2 



(3) 



which corresponds to r(q) = — ln[a q + ft 9 ]/ In 2 and can be obtained from a 
generalized binomial multifractal model [13], see also [4,8]. The values of the 
two parameters a and b are also reported in the figures. The results for all 
rivers can be fitted surprisingly well with only these two parameters (see Fig. 
3). Instead of choosing a and b, we could also choose the Hurst exponent h(l) 
and the persistence exponent h(2). From knowledge of two moments, all the 
other moments follow. 

This surprising result does not hold for the precipitation records. As can be 
seen in Figs. 4(c) and 4(f) there are stations where Eq. (3) cannot describe 
the multifractal scaling behaviour reasonably well. According to Rybski et al., 
Eq. (3) is appropriate only for about 50 percent of the precipitation records 
[14]. 
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In the generalized binomial multifractal model, the strength of multifractal- 
ity is described by the difference of the asymptotical values of h(q), Aa = 
h(—oo) — h(oo) = (ln6 — In a)/ In 2. We note that this parameter is identical 
to the width of the singularity spectrum /(a) at / = 0. Studying 41 river 
runoff records [13], we have obtained an average Aa = 0.49 ± 0.16, which 
indicates rather strong multifractality on the long time scales considered. For 
the precipitation records, on the other hand, the multifractality is weaker. The 
average is Aa = 0.29 ± 0.14 for 83 records [14]. 



Our results for h(q) may be compared with the different ansatz h(q) = 1+H' — 
Ci(g a ' _1 — l)/(a' — 1) with the three parameters H', C-y, and a' (LS ansatz), 
that has been used by Lovejoy, Schertzer, and coworkers [6,7] successfully to 
describe the multifractal behaviour of rainfall and runoff records for q > 0. 
A quantitative comparison between both methods is inhibited, since here we 
considered only long time scales and used detrending methods. We like to note 
that formula (3) for h(q) is not only valid for positive q values, but also for 
negative q values. We find it remarkable, that for the runoff records only two 
parameters were needed to fit the data. For the precipitation data, one needs 
either three parameters like in the LS ansatz or different schemes. 



In summary, we have analyzed long river discharge records and long precipita- 
tion records using the multifractal detrended fluctuation analysis (MF-DFA) 
and the wavelet transform modulus maxima (WTMM) method. We obtained 
agreement within the error margins and found that the runoff records are 
characterized by stronger multifractality than the precipitation records. Sur- 
prisingly, the type of multifractality occurring in all runoff records is consistent 
with a modified version of the binomial multifractal model, which supports 
the idea of a 'universal' multifractal behaviour of river runoffs suggested by 
Lovejoy and Schertzer. In contrast, according to [14], several precipitation 
records seem to require a different description or a three-parameter fit like the 
LS ansatz. The multifractal exponents can be regarded as 'fingerprints' for 
each station. Furthermore, a multifractal generator based on the modified bi- 
nomial multifractal model can be used to generate surrogate data with specific 
properties for each runoff record and for some of the precipitation records. 
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